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We examine single-frequency optical schemes for species-selective trapping of ultracold alkali- 
metal atoms. Independently addressing the elements of a binary mixture enables the creation of an 
optical lattice for one atomic species with little or no effect on the other. We analyze a "tune-in" 
scheme, using near-resonant detuning to create a stronger potential for one specific element. A 
"tune-out" scheme is also developed, in which the trapping wavelength is chosen to lie between two 
strong transitions of an alkali-metal atom such that the induced dipole moment is zero for that 
species but is nonzero for any other. We compare these schemes by examining the trap depths and 
heating rates associated with both. We find that the tune-in scheme is preferable for Li-Na, Li-K, 
and K-Na mixtures, while the tune-out scheme is preferable for Li-Cs, K-Rb, Rb-Cs, K-Cs and 
39 K- 40 K mixtures. Several applications of species-selective optical lattices are explored, including 
the creation of a lattice for a single species in the presence of a phononlike background, the tuning 
of relative effective mass, and the isothermal increase of phase space density. 

PACS numbers: 33.80.Ps, 03.75.Nt, 03.75.Ss, 39.25.+k 
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I. INTRODUCTION 



As the field of ultracold atoms research enters its ado- 
lescence, experiments are increasingly including more 
than one element or isotope. Dual-species experi- 
ments offer possibilities for creating heteronuclear polar 
molecules [l[ , sympathetic cooling & S S IE| , and in- 
vestigating Bose-Fermi mixtures [3, |j, la, |6( which may 
provide opportunities for studying boson-mediated su- 
perfiuid states @, [1] ■ More recently, experiments involv- 
ing up to three atomic species have been implemented 
for sympathetic cooling of two fermionic species [l(| ■ 

Dually degenerate experiments have so far used exter- 
nal trapping potentials common to both atomic species. 
A species-specific trapping potential would add a degree 
of freedom to improve sympathetic cooling |ll| . tune ef- 
fective mass, or create a lattice for one species in the pres- 
ence of a background reservoir. Though careful selection 
of internal atomic states can provide differential mag- 
netic trapping, optical far-off resonant traps (FORTs) 
and magnetostatic traps are not species-specific. 

Species-selective adiabatic potentials have been pro- 
posed [HI] and demonstrated [IH in the case of 87 Rb- 
40 K, where the Lande factors |<7f| are distinct. A radio- 
frequency transverse field can be resonant with only one 
of the two species, selectively deforming its dressed po- 
tential. Onofrio and co-workers [ll|> [l4| propose using 
two overlapping FORTs at frequency detunings far above 
and below the dominant ground state transitions of both 
species in a two-species mixture. The confinement of 
each species can be chosen independently by individu- 
ally adjusting the intensity of the two beams used to 
create the trap. Unfortunately, neither of these schemes 
lends itself to a uniform lattice potential for atoms: the 
radio-frequency scheme fails because it is limited to one- 
dimensional periodic potentials and the two-frequency 
balancing because lattice periodicity depends on wave- 
length 



In this paper we discuss the generation and application 
of species-specific optical lattice potentials. Our motiva- 
tion is the strong analogy between atoms in optical lat- 
tices and electrons in crystalline solids. Cold bosons in 
lattice potentials can be used to explore strongly interact- 
ing many-bod y p hysics, such as the supcrfluid-insulator 
transition [lq.ll7|. At sufficiently low temperatures, cold 
fermions in lattices [H, Ojj] might be able to address open 
questions about the ground state of the Hubbard Model 
20]. Unlike crystal lattices, optical lattices do not sup- 
port the lattice vibrations responsible for many physical 
phenomena. However, phonon mediation between neu- 
tral atomic fermions could arise in the presence of a con- 
densed bosonic s pec ies capable of sustaining phonon-likc 
excitations |7| la, |2lJ, |22j . Although Bose-Fermi mixtures 
were recently loaded into optical lattices [H, [24j], the 
lattices confined bosons as well as fermions. A lattice- 
induced increase in the effective mass of the bosons re- 
duces the speed of sound in the condensate, and thus the 
mediating effects of the phonons. If the fermions were 
confined in a lattice without effect on the bosonic cloud, 
this mediation would be less inhibited [251 ]. 

We discuss two approaches to species-specific optical 
potentials, both of which involve only a single frequency 
of laser light. The first approach is to tune the laser 
wavelength close to the atomic resonance of one species, 
making its induced dipole moment much stronger than 
that of any other atomic species present. We refer to 
this strategy as the "tune-in" (TI) scheme. A second ap- 
proach exists for atoms, such as alkali metals, with an 
excited state fine structure splitting. Between the reso- 
nances of the doublet, a wavelength can be chosen such 
that the induced dipole moment is strictly zero. We refer 
to this strategy as the "tune-out" (TO) scheme. Both ap- 
proaches allow for the creation of a species-specific optical 
lattice with a tunable relative potential strength between 
species. 

In the following sections we consider the relative mer- 
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its of the TI and TO schemes. We focus our attention 
on Bose- Fermi combinations throughout, pa ying special 
attention to mixtures including 87 Rb O Ea E2 EE El, 
l30l . [3~0 | . Section [TT] presents the methods used to calcu- 
late potential energies and heating rates of atoms in a 
light field. We present the TI and TO approaches in 
Sec. QUI and quantify their relative merits. For several of 
the stable alkali-metal isotopes, we calculate the tune-out 
wavelengths, scattering and associated heating rates, and 
the trapping potentials imposed upon a second species 
of alkali-metal atom. Section IIVI discusses the role of 
interactions, both in limiting the distinct addressability 
of one species and in allowing thermalization between 
two species. In Sec. [V] we discuss applications of the 
species-specific dipole potentials for cooling, engineering 
mobility, and adding phonons to optical lattices, before 
concluding in Sec. I VII 



II. METHODS 

We consider the light-atom interaction in the limit of 
a small excited-state fraction. An electromagnetic field 
induces an electric dipole potential on a neutral atom, 
calculated in Sec. Ill Al which can be used as a trapping 
potential for ultracold atomic samples. We consider the 
residual effect of spontaneous emission in mi Bl 



A. Dipole potential 

An atom in a ground state \g) will experience a po- 
tential shift due to coupling by the light field to the ex- 
cited states |e). We calculate the sum of these shifts on 
each state \g} 7 including the counter- rotating term, using 
second-order perturbation theory: 
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where u>l is the laser frequency, Tiu) eg is the energy dif- 
ference between \e) and \g), d is the dipole operator, e is 
the polarization of the light, and I is the light intensity. 
In the case of atoms in a weak magnetic field, we use the 
matrix elements defined in the \F, mp) basis, where F 
is the total angular momentum and mp is the magnetic 
quantum number. 

We will consider only alkali-metal atoms, which have 
two dominant ns — > np transitions due to the fine struc- 
ture splitting. Using nomenclature established by Fraun- 
hofer for the 3 2 S 1/2 -> 3 2 P 1/2 and 3 2 S 1/2 -> 3 2 P 3/2 
transitions in sodium, we label the corresponding lines in 
each of the alkali metals D\ and D 2 , respectively. Spin- 
orbit coupling splits each excited state by a frequency 
Aps = ^d 2 ~ ^Di, while each ground and excited state 
is further split by the hyperfine interaction Ahfs an d 
A^ps, respectively. The atomic data used for Eq. ((T|) 
are the measured linewidths and line centers of the D\ 



and D 2 li nes, and the ground and excited state hyperfine 
splittings [321 ]. 

Transitions to higher excited states ns — > (n + l)p are 
neglected by our treatment. When detuned within Aps 
of the ns — > np transition, the relative magnitude of the 
ns — > (n + l)p shift is less than 2 x 1CP 5 for Cs and 
7 x 1(T 8 for Li. 

As Eq. |T]) requires a sum over several states and 
knowledge of individual matrix elements, it is useful to 
have an approximate but simpler expression for U g . If 



the detunings A e 



<jJl 



are small compared to 



ApSi but large compared to the excited state hyperfine 
splittin g A ^^q , an approximate expression for the dipole 
shift is I33| 
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where P = 0, ±1 for w, polarization, respectively, gp 
is the Lande factor, A^) is the detuning from the D]/ 2 ) 
line, ujq = (lud 1 + 2o;d 2 )/3 is the line centre weighted by 
line strength, and T = (Td 1 +Td 2 )/2 is the average of the 
Di and D 2 linewidths. Since Tjj 2 /Td 1 »1 + 3Aps/^o, 
one can expect an accuracy between ±7% for Cs and 
±0.003% for Li. At small detunings, we empirically find 
that Eq. deviates from Eq. ([T]) by >1% for detunings 
min{|Ai|, |A 2 |}/2tt < 1.5-/M GHz, where M is the mass 
in atomic units. 

As an approximate form, Eq. @ neglects the counter- 
rotating terms. For |A| < Aps, the strength of the 
counter-rotating contribution is at most Afs/2wl rela- 
tive to the contribution of one near-resonant dipole tran- 
sition. Thus the neglected shift is at most -2% for Cs, 
and -0.001% for Li. 



B. Heating Rate 

Detuning and intensity of optical traps must be chosen 
with consideration of the incoherently scattered trapping 
light that heats the atoms. For each state \g), we quantify 
the rate of scattering in the low saturation limit, 
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where T e is the natural linewidth of the j -t e transi- 
tion, and A eg = lul — Lu eg . The rate of scattering of 
photons can be converted to an average heating rate, 
H g = ^Eji'jg, where E R = 7i 2 fc 2 /2m is the recoil energy, 
k = lul/c, and c is the speed of light. 



III. SPECIES-SPECIFIC DIPOLE TRAPPING 
SCHEMES 

It is not surprising that optical traps can be species- 
specific, given that optical resonances are unique to 
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Target 


P 


7 Li 


23 Na 


Spectator 

39 K 


87 Rb 


133 Cs 


6 Li 


oo 


0.00134 


7.77xl0 -4 


-0.0381 


-1.20 


-8.45 




100 


2.66xl0~ 5 


0.281 


0.220 


0.239 


0.347 




10 


3.92xl0~ 6 


2.54 


2.30 


2.49 


3.57 


40 K 


oo 


4.28xl0~ 7 


5.77xl0 -4 


0.188 


-9.03 


-25.8 




100 


3.05 


5.78 


3.64xl0 -5 


0.251 


1.39 




10 


29.2 


53.0 


4.65xl0 -4 


3.27 


18.8 



TABLE I: Sustainability s of two-species mixtures for tune-out and tune-in schemes, in units of seconds, using Eq. ([5]). The 
second column indicates the selectivity. Rows with p — oo are calculated using the TO scheme while rows with p = 100 and 
p = 10 use the TI scheme. For the lighter spectators, the TI scheme has higher |s|. The heaviest elements and isotope mixtures 
favor the TO scheme. 



atomic elements and isotopes. However, most species- 
specific optical traps, such as magneto-optical traps, are 
tuned to within a few linewidths of resonance, which is in- 
compatible with quantum degenerate ensembles. At low 
temperatures and high density, any gain in trap depth 
close to resonance must be balanced against the heating 
due to unwanted light scattering (Sec. Ill Bp . 

In the subsections below we consider two-species mix- 
tures. The goal is to apply a dipole force to the "target" 
species while inducing as little potential as possible on 
the second species, which we will call the "spectator". 
We define the "selectivity" as 



Ut 

Us 



(4) 



where U t r s -\ is the potential induced on the target (spec- 
tator). 

As discussed in Sec. UH in the low saturation limit, 
both the induced dipole potential and the heating rate 
are proportional to intensity. We define the intensity- 
independent ratio 



s = 



(5) 



to be the "sustainability" , where -Ht(» is the heating rate 
of the target (spectator) . The absolute value of s sets the 
scale for possible trapping time [33] • The laser frequency 
will be chosen to maximize both p and |s|. 



A. The tune-in scheme 

The simplest selective potential is one in which the 
laser is tuned close to a resonance of the target (see 
Fig. [I}. Since scattering rate is inversely proportional 
to the square of the detuning, we consider only heating 
of the target such that sti — > U t /H t . Considering Eq. ((2|) 
in the limit |A t | <C |A S |, a simple estimate is 



Sti 



3h 



(6) 



The choice of detuning A t will depend upon the desired 
selectivity. Assuming that the spectator is far-detuned, 
A s ps (u> ,t - Wo,a) and 



Pti 



2r t wX, N,t 



3r, 




(7) 



Together, Eqs. ([6]) and (0 explicitly show the opposing 
dependence on A t and the necessary trade-off between 
selectivity and sustainability. 

Table U shows the sustainability for bosonic species as 
spectators and fermionic species as targets, calculated 
using Eqs. ([I]) and ((3]) for tune-in selectivities of 100 and 
10. The inverse scaling predicted by Eqs. © and ([7]) 
is observed: Sti drops by approximately a factor of ten 
when pti is increased by a factor of ten. The product 
of Sti and Pti in Table |T] varies between 22 s and 36 
s for 6 Li mixtures and between 140 s and 530 s for 40 K 
mixtures, whereas Eqs. ([6]) and ||7J) would predict a ranges 
of 8 s to 21 s and 64 s to 132 s, respectively, excluding 
isotope mixtures from the comparison. Finally, we note 
that 23 Na- 40 K is the optimal tune-in mixture. 



B. The tune-out scheme 

The tune-out wavelength scheme exploits the charac- 
teristic doublet structure of the alkali-metal atoms. By 
choosing a wavelength that lies between the two strongest 
transitions, the large negative energy shift of the D2 line 
is balanced against the large positive energy shift of the 
Di line (see Fig. [T|). This atom becomes a spectator, 
while any other species feels the shift induced by the 
laser and becomes a target. Since the potential shift on 
the spectator can be zero, the selectivity of the tune-out 
approach is infinite. 

The laser frequency, wtoi a t which U s = 0, is de- 
termined numerically using Eq. |T]) with all \F, TOf) ex- 
cited states in the D\ and D2 manifolds. Table [TT] shows 
the tune-out wavelength for all 87 Rb ground states with 
7r, linear [35| . <r + , and a~ polarizations. An approxi- 
mate expression for this wavelength can be derived from 
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' 775 780 785 790 795 800 



laser wavelength (nm) 

FIG. 1: Energy shift as a function of wavelength for 87 Rb 
in the \F, tuf) = 1 2, 2) state, under linear polarization, for 1 
mW/cm 2 . This general structure will arise for each of the 
alkali-metal elements, with the divergences located at the Di 
and Z?2 lines. The tune-in scheme for a 87 Rb target is in- 
dicated by the square marker on the blue-detuned branch of 
the potential energy curve. The round marker indicates the 
position of the tune-out wavelength, where the energy shift is 
zero. Here, 87 Rb is the spectator. 

Eq. (0), giving 

1 + Pg F m F 

"TO — ^0 g ^FS- \°) 

For 87 Rb, Eq. (JSJ) predicts tune-out wavelengths of 785.10 
nm, 787.54 nm, 790.01 nm, and 792.49 nm for Pgpmp = 
-1, -0.5, 0, and 0.5, respectively. Comparing these values 
to the results of Table HT1 we see the that approximations 
are accurate to 0.04 nm or better. Equation §E§ also 
gives wto = for the case of PgFmp = 1, which 
is inconsistent with the assumption that |Ai| 3> Ahfs- 
In fact, since PgFmp — 1 corresponds to a dark state 
with respect to the Di excitation, there is no tune-out 
wavelength for this case. 

We note that in Table [TT] the tune-out wavelength for 
linear polarization is nearly independent of the choice 
of ground state, in contrast to er + or o~ polarizations. 
Given this independence, we calculate the tune-out wave- 
lengths using Eq. ((T|) for several of the common alkali- 
metal isotopes in their stretched ground states under lin- 
ear polarization (Table ITTT|) [361 ]. 

Tables UU and IIIII also give the scattering rates per 
unit intensity (scattering cross-section) for the specta- 
tor species at the tune-out wavelength. Due to the large 
dispersion of fine splitting among the alkali metals, scat- 
tering cross-sections vary by over six orders of magnitude. 
To understand the implications for trapping, we need to 
consider the sustainability, s, of various pairs of atomic 
species. Table U shows sto i n the same mixtures for which 
sti is calculated. Spectators with larger scattering cross- 
sections at Ato have lower sustainability. In addition 
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9.1xl0~ 6 
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TABLE II: Tune-out wavelengths and scattering rates in 
87 Rb, for select polarizations and all ground states. Two 
states have no tune-out wavelength, because the D\ line has 
no F=3 excited state. 



Element \F, tuf) Ato (nm) Jsc/I 







Eq.[T] 


Eq.[2] 


(cm 2 /mj) 


6 Li 


|3 3\ 
1 2' 21 


670.99 


670.99 


2.8 


7 Li 


|2,2) 


670.97 


670.97 


2.4 


23 Na 


|2,2) 


589.56 


589.56 


2.0xl0" 3 


39 K 


|2,2) 


768.95 


768.95 


1.4xl0~ 4 


40 K 


1 9 9\ 
13' 2/ 


768.80 


768.80 


1.7xl0~ 4 


87 Rb 


|2,2) 


790.04 


790.01 


9.1xl0" 6 


133 Cs 


|4,4) 


880.29 


880.06 


1.5xl0~ 6 



TABLE III: Tune-out wavelengths and scattering rates for 
various elements. We have assumed linear polarization and 
stretched states and show a comparison of Eqs. |T} and ([2]). 
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to these data, we note that the 133 Cs- Rb spectator- 
target combination gives the highest possible sustainabil- 
ity among alkali metals in the tune-out scheme: s = 34 s 
at a tune-out wavelength of Ato = 880.29 nm (not shown 
in tables). 

As an example, consider making a lattice for 40 K only, 
leaving 87 Rb unaffected by the lattice and confined only 
by a background magnetic trap or FORT. If we con- 
sider the mixture in a three-dimensional lattice of ar- 
bitrary depth, the 40 K potential shift is 1.54 x 10~ 5 fiK x 
[/(mW/cm 2 )]; with beams of 100/jm waist, the potential 
shift is 98 nK/mW. If we require the target to experi- 
ence a lattice that is 8Er deep [I?}, where Er is the 
recoil energy, we find a 87 Rb heating rate of 98 nK/s. 

C. Discussion: tune-in vs. tune-out 

With the data of Table U we can evaluate the practi- 
cality of both the tune-in and tune-out schemes for Bosc- 
Fermi mixtures of neutral alkali metal atoms. Since the 
scattering rate in the tune-out scheme can be smaller 
for elements with larger fine structure splittings, this ap- 
proach is better suited to more massive elements [38]. 
In particular, the tune-in scheme is preferable for Li-Na, 
Li-K, and K-Na mixtures, and for applications requir- 
ing selectivity of less than 10:1. The tune-out scheme is 
preferable for Li-Cs, K-Rb, and K-Cs mixtures when the 
selectivity required is greater than 10:1, and for Li-Rb 
mixtures at selectivity of greater than 20:1. 

An isotope mixture of potassium could be compatible 
with the tune-out scheme, where s > 100 ms. Isotope- 
specific manipulation within a lithium mixture is less 
practical due to sustainabilities of 1 ms or less. 

Other factors may also influence whether the tune-in 
or tune-out approach is preferred. For experiments with 
time scales that are slow compared to the thermalization 
rate (discussed further in Sec. IIVB[) . it may be prefer- 
able to heat the minority species and allow for sympa- 
thetic cooling. If the reservoir of spectator atoms is large, 
the tune-in scheme might be preferred since extra energy 
due to near-resonant heating would be transferred to the 
reservoir. For experiments on time scales fast compared 
to the thermalization rate, the tune-out scheme might be 
preferred even if s is smaller, since spectator heating will 
not affect the target. 

Finally, we note that if a third species is included, it 
will play a spectator role in the tune-in scheme and a 
target role in the tune-out scheme. 



accordingly. A more relevant characterization of trapping 
and dissipation is 

C/ pk 

Sglobal = (H t )l(H s y (9) 

where the brackets indicate expectation values and [7 t max 
is the maximum potential depth of the target. In the 
case where only one heating rate dominates, we define a 
quantity [3 to represent the enhancement of sustainability 
due to the lattice: 

/Jto.ti^^-t^ , (10) 

where the expectation value is taken for the spectator 
species in the tune-out case, and for the target in the 
tune-in case. In a one-dimensional lattice, the intensity 
is I(x) — Iosm 2 (kx), where k — 2tt/\l. 

For the tune-out scheme, neither the position nor the 
motion of the spectator is coupled to the standing wave 
intensity. We can assume that all intensities are equally 
sampled and use the average intensity; (I) s t = I /2 to 
give /3 T o = 2. 

For the tune-in scheme, the target is the species for 
which the heating rate is dominant and for which the 
lattice enhancement of the sustainability must be taken 
into account. We first consider (3 in the case of a local- 
ized atom. Classically, we expect deeply bound states 
of a fe/ue-detuned lattice to avoid the intensity maxima 
where incoherent scattering is high, such that the inten- 
sity seen by the atom is less than the average. Heating is 
therefore reduced and /3ti is increased. For an atom local- 
ized to a single site by interactions, we can use Wannier 
states to calculate the expectation value of the intensity 
(see Fig. [2]). For lattice depths 17 = U? k /E R > 1, these 
states are approximately harmonic oscillator states with 
oscillator frequencies Who = 2V 'EuUt/fi. In this case, 
/?ti ~ 2y / r/; in an arbitrary number of dimensions, d, the 
advantage is increased as /3ti ~ 2 d r\ d l' 1 for rj 3> 1. 

Considering that classical oscillators spend more time 
at turning points than at potential minima, blue tune-in 
lattices may not always have (3to > Pti- For a delocal- 
ized q = 1.5hk Bloch function (in the first excited band), 
where q is the quasi-momentum, we show f3 versus 77 in 
Fig. [2] For 77 < 10, confinement effects give /3ti < 2. 

IV. INTERACTIONS BETWEEN SPECIES 
A. Mean field interactions 



D. Lattice confinement effects 

In Sees. A-C we have discussed the ratio of induced 
potential to heating rate as a local quantity. However, 
the spatial dependence of intensity in an optical lattice 
means that moving atoms may not sample all laser inten- 
sities uniformly and the heating rate should be modified 



Interactions between elements may couple the envi- 
ronment felt by the target to the spectator, spoiling the 
species-specific potential shaping of both schemes. If tar- 
get atoms are trapped in a lattice in the presence of a 
background spectator species, a periodic interaction po- 
tential could arise for the spectator due to its interac- 
tion with the periodically modulated density of the tar- 
get. A mean- field approach is used to make a simple 
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5 10 15 

Lattice Depth U^/E K 



20 



FIG. 2: A standing wave of laser light increases the ratio of 
confinement depth to scattering rate, modifying the sustain- 
ability s calculated for uniform intensity. The enhancement 
fori is shown as a function of r] (trap depth) for a localized 
state (solid line), a q = Bloch state in the lowest band, 
and a q = 1.5hk Bloch state in the first excited band (both 
dashed). At high lattice depths, both wavefunctions in the 
lowest band approach the tight binding limit (dotted line). 
For comparison, foro = 2 for all lattice depths. 



estimate of the magnitude of this effect. The interac- 
tion potential of the spectator due to the target is given 
by Uint.s — (4:Trh 2 a s t/ n s t)nt, where a st is the scattering 
length between species, n t is the density of the target, 
and n s t = 2m s mt/(m s + mj) is the reduced mass. As 
Ut/Ui n t.s — > 1, interaction effects begin to degrade selec- 
tivity. Estimating the three-dimensional density of the 
target at the centre of a single lattice site using the har- 
monic oscillator ground state with frequency Who and as- 
suming one atom per site, we find an interaction-limited 
selectivity 



1/4 



327r|a s t|m 3 / 2 E 1 ^ 2 



(11) 



In the case of the 6 Li- 87 Rb target-spectator mixture, 
the interaction- limited selectivity is p max = 3. lry 1 / 4 [3l| . 
while for the 40 K- 87 Rb target-spectator mixture, it is 
Pmax = O.2277 1 / 4 . To be in a regime where interaction 
effects might be ignored, we require p max 3> 1 which 
gives r) > 0.01 for 6 Li- 87 Rb and 77 > 400 for 40 K- 
87 Rb. While for the lithium target, this condition is quite 
reasonable, the large lattice depths required for potas- 
sium would completely localize the atoms to individual 
sites and prevent the exploration of interesting tunneling- 
driven physics. Overcoming this interaction limitation 
and achieving the high selectivity discussed in Sec. IIIII 
may be possible by tuning a magnetic field to a value 
where aRt>K ~ near a Feshbach resonance (39j . For in- 
stance, with depth rj = 8 and p max > 10, the scattering 
length between species is limited to |a st | < 7.5ao, where 
cto is the Bohr radius. 



Where species selectivity is not the goal, this 
interaction-induced periodic potential for a second 
species could be used alternative lattice potentials. Such 
potentials are non-sinusoidal, do not involve a Stark shift, 
and may have a dynamic structure if target atoms are 
mobile. The strength of this induced potential could be 
controlled through the interaction strength, such as at a 
Feshbach resonance, as described above. 



B. Inter-species thermalization 

An understanding of the thermalization between 
species is important when considering heating or cool- 
ing in the trap. The rate at which energy is transferred 
between the species will be relevant in setting the time 
scales on which adiabatic experiments can take place. 
Thermometry is also possible if there is good thermal- 
ization between species; a high density target could be 
confined within a species-specific lattice while the specta- 
tor remains extremely dilute and thus at lower quantum 
degeneracy, where temperature is more easily measured. 
Conversely, thermal isolation could be useful in shielding 
one species from the spontaneous heating in the other. 

In the classical limit, the thermalization rate is propor- 
tional to the collision rate of the atoms in the trap, given 
by jcoii — nav, where n is the overlap density, a the scat- 
tering cross-section, and v the relative velocity between 
species. Random collisions act to equilibrate the system 
and the rate of rethermalization is jtherm ~ 7co/;/2.7 [40j . 

For a degenerate mixture of bosons and fermions, the 
classical picture of scattering breaks down and the rate 
of thermalization decreases as the fermionic system be- 
comes more degenerate. An estimate of the sympathetic 
cooling of a uniform system in the degenerate regime us- 
ing the quantum Boltzmann equations gives the rate of 
change in the degeneracy (4l| : 



d_ (T_ 

It V Tf 



6C(3) 



~fcoll 



T 



(12) 



where j' co n = (3/8)nB<JVF is the collision rate between 
species, with the scattering cross-section a — Ai:a 2 BF , 
Fermi velocity, v F = fi,(67m F ) 1/3 /m F , C(3) ~ 1.20206 is 
the Riemann zeta function and Tp is the Fermi temper- 
ature. Though other assumptions [42j of this treatment 
are not valid for the systems we consider in this paper, we 
use this expression to determine an order of magnitude 
for the thermalization rate. 

Taking boson density tib — 1 x 10 14 cm~ 3 , fermion 
density np = 1 x 10 13 cm -3 , and using the mass for 
40 K, we calculate that d/dt(T/T F ) w -350s" 1 (T /T F ) 2 , 
which gives a temperature relaxation time of « 30 ms at 
T = O.lTp. This is an order of magnitude larger than the 
the classical expectation for the rethermalization time of 
w 2 ms for particles moving at the Fermi velocity. Thus 
for rapid experiments in the deeply degenerate regime, 
thermal contact is essentially broken, allowing, for in- 
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stance, the target to be unaffected by the heating of the 
reservoir in the tune-out scheme. 



V. APPLICATIONS OF SPECIES SPECIFIC 
TRAPPING 

A. Isothermal phase space increase 

Two-species mixtures can be used to realize various 
cooling schemes. For example, dark state cooling by su- 
pcrfiuid immersion is discussed in Ref. (43|. We present 
two simple cooling scenarios in which the presence of an 
uncompressed spectator allows the target species to be 
compressed with negligible temperature increase but sig- 
nificant improvement in phase space density. In both 
cases, target atoms are first compressed isothermally, 
spectator atoms are then removed from the trap, and 
finally, the target atoms are decompressed adiabatically. 

We consider a species-selective single well dipole trap 
and a species-selective lattice. A closed cycle is described 
for both scenarios with atoms beginning and ending in a 
FORT. We assume that the heat capacity ratio between 
spectator and target is infinitely large and that thermal- 
ization in the spectator and between the spectator and 
target species is faster than any other time scale consid- 
ered. The latter assumption may restrict fermion cooling 
to the non-degenerate regime (see Sec. IIVBI) . 

In the first scenario, schematically represented in Fig. 
[31 a single species-specific beam crosses a FORT. Adia- 
batic cooling reduces temperature in proportion to the 
ratio of average trapping frequencies, that is, Tj/Tf — 
uji/ojf, where i,f indicate "initial" and "final". As an 
example, consider 40 K- 87 Rb and 6 Li- 87 Rb mixtures, con- 
fined by a 1064 nm, 500 mW single-beam FORT with a 
1/e 2 radius of 20 fim and a corresponding trap-averaged 
harmonic oscillator frequency u>i = 2ir x 540 Hz for K 
and Wi = 2tt x 950 Hz for 6 Li. A 500 mW, 50 fim waist 
beam at Ato (see Table llllj) is turned on perpendicular 
to the FORT to compress the fermions in a trap with 
frequency tJ f = 2tt x 3040 Hz for 40 K and uJ f = 2tt x 3810 
Hz for lithium. After providing sympathetic cooling dur- 
ing the compression of the target atoms, the rubidium is 
ejected from the trap by temporarily removing the FORT 
or by using a resonant pulse of light. The species-specific 
trap is then adiabatically ramped down and turned off, 
leaving the fermions in the FORT at a temperature ap- 
proximately 5.7 (4.0) times colder than when they began 
for potassium (lithium). Though this is a modest change 
in temperature, the phase space density in a harmonic 
trap is proportional to the inverse cube of the temper- 
ature, indicating a factor of 180 increase in phase space 
density for potassium and a factor of 60 for lithium. 

In the second scenario, we consider a 3D lattice created 
by the tune-out wavelength 0|. The target-specific lat- 
tice is ramped on until peak lattice intensity is reached. 
The spectator is evaporatively cooled and ejected by re- 
ducing the spectator trap depth. The lattice is then 




FIG. 3: A cooling procedure using species-specific trapping; 
one- dimensional trap shape represents three-dimensional 
trap-averaged shape. The solid line represents a FORT, the 
dashed line the species-specific trap; open circles are the spec- 
tator species and closed circles are the target, (a) Two species 
are trapped in a FORT; (b) the species-selective beam is 
turned on, compressing and heating the target species; (c) the 
target species rethermalizes with the spectator; (d) the spec- 
tator is removed; (e) the target is adiabatically decompressed 
to a lower temperature and transferred to a FORT. 



ramped back down isentropically, leaving the target in 
the initial trap with an entropy and temperature limited 
by the "plateau entropy" discussed in [45|. As shown 
there, a target of fermions with unity filling has an en- 
tropy plateau of zero, which would suggest no lower limit 
to the achievable temperature. 

An important limitation of these schemes will be the 
competition between adiabaticity and heating. For an 
ideal gas, the condition of adiabaticity requires any 
changes to take place in times longer than the inverse 
of the smallest trap frequency. Using the numbers given 
in the 40 K- 87 Rb example of the crossed dipole trap in 
the tune-out scheme, we find that adiabaticity requires 
a relaxation time of 35 ms, at an intensity of 3.5 xlO 7 
mW/cm 2 , yielding a heating of 1.3 /iK during this time, 
which sets a lower bound for the temperature attainable 
in this scheme. Another possible limitation of the cooling 
schemes is the interspecies thermalization, which limits 
the speed of the isothermal step (see discussion in EUVBI) . 



B. Phonons 

Unlike crystal lattices of solids, optical lattices do not 
support phonons. Since these quasi-particles play a lead- 
ing role in the physics of condensed matter, it is of in- 
terest to introduce phononlike excitations into a system 
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of ultracold atoms in an optical lattice. The boson- 
mediated interaction between fermions has been stud- 
ied in both uniform 0, H SH, HI and lattice systems 
[25l l46l l47j . Little theoretical work has been done on a 
system in which only the fermions are confined to the lat- 
tice [H)]. Here, the background bosons are free to interact 
with the fermions and with one another. If the bosons 
are degenerate, the condensate can sustain phonon-like 
excitations and allow for boson-mediated interactions be- 
tween fermions on spatially separated lattice sites. 

For the phonons in the condensate to play a role in 
mediating interactions, their spatial extent must exceed 
the lattice spacing. The healing length of a uniform bose 
condensate, £ = (8na) -1 / 2 where n is the condensate 
density and a is the scattering length between species, 
sets the relevant length scale. Species-specific trapping 
allows the superfluid bosons to remain at low density 
while fermions are tightly bound in the lattice, thereby 
maximizing the range of the mediated interaction. A 
finite selectivity does not prevent mediation of interac- 
tions, since the bosons remain superfluid at depths less 
than the Mott-insulator superfluid transition [17], per- 
mitting both tune-in and tune-out schemes to be used 
for this application. 




10 20 30 



Lattice Depth (J7 t pk /£ R , t ) 



C. Effective mass tuning 

An optical lattice can be used to change the effective 
mass of the atoms in it, allowing for the tuning of ex- 
perimental parameters including interaction strength [49[ 
and tunneling rate (50| . which can be used, for example, 
to explore different regimes of collective dynamics [5l| . 
The effective mass of a wave packet centered at quasi- 
momentum q c is 



m* (q c 



d 2 E 



(13) 



where E(q) is the band energy. Figure Ufa) shows 
the effective mass for q c = and q c = 0.5hk in 
a one-dimensional optical lattice potential U (x) = 
U pk sin 2 (kx). 

For deep lattices, an approximate form for the q c = 
case is m* = h 2 k 2 /2ir 2 J, where J is the tunnelling energy 
giving an effective mass enhancement 



to, 47T 3 / 2 ?? 3 / 4 



(14) 



The ratio of effective masses for the target and spectator 
can be estimated from Eq. (TJJ 



1/4 



mz 



l t „ 6XP { 2 ^* V ~ V m */P m t) } ( mt n 
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Figure BJb) shows the effective mass ratio for the case 
of a 6 Li target and a 133 Cs spectator. In particular, it 



FIG. 4: (a) The ratio of the effective mass to bare mass of 
the target species is shown as a function of trap depth r/t, for 
both q — and q c = 0.5hk, as labeled. The tight binding 
approximation (dashed line) approaches the exact calculation 
for rjt 3> 1. (b) The effective mass ratio is shown for the case 
of a 6 Li target and a 133 Cs spectator. Whether the ratio in- 
crease or decreases depends on the selectivity p of the lattice. 
The critical point is when p — m s /m t (= 22.2 in this case), 
as explained in the text. The tight binding approximation of 
the effective mass ratio is also shown (dashed line). 



is striking that for p — 10, the ratio decreases with lat- 
tice depth, while for p — 100, the ratio increases with 
lattice depth. The critical selectivity is well predicted 
by Eq. (fT5|) ; at p = m s /m t , the effective mass ratio is 
independent of lattice depth. 

Both the tune-out and tune-in schemes provide the 
means of choosing selectivity. Several tune-in selectiv- 
ities are shown in Table HJ the tune-out selectivity can be 
chosen simply by choosing a wavelength slightly different 
from Ato- In the example used here, the sustainability 
for p = 10 is better for the tune-in scheme, but tuning to 
to*^ > m^ s at moderate lattice depths requires a p — 100, 
for which the tune-out scheme has higher sustainability. 



VI. CONCLUSIONS 

We have discussed how the choice of wavelength used 
to create an optical lattice can tune its selectivity be- 
tween elements or isotopes. This control should increase 
the range of parameters that can be explored in multi- 
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species ultracold atom experiments. The tune-out wave- 
length scheme allows for the complete cancellation of the 
trapping potential for one species while providing a con- 
fining or lattice potential for any other species in the sys- 
tem. This scheme will work best using the heavier alkali- 
metal atoms, Rb and Cs, as the spectator elements, and is 
most successful for the 40 K- 133 Cs fermion-boson mixture. 
The alternative tune-in scheme uses a near-detuned opti- 
cal potential, creating a much stronger potential for one 
element than the others, without the ability to strictly 
cancel the potential for one element. Mixtures involving 
the Li, Na, and K as spectators are most compatible with 
this approach, where the 40 K- 23 Na is the most promising 
fermion-boson mixture. 

The power of species selection is in its use to engi- 
neer specific lattice environments for the atoms, includ- 
ing adding a bosonic background to mimic the phonons 
present in solids and tuning the relative effective mass 
of two species. Applications for which experiments are 
rapid compared to the sustainabilities calculated in Sec. 



IIIII are especially promising. Species selection enables 
cooling in a two-species mixture, and in the case of 
fermions trapped in a lattice, reduces the temperature 
of fermions as they are loaded into a lattice, in contrast 
to current experimental realizations [45|, l53| . 
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